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Abstract. In the joint work [RZ] of T. Rivoal and the author, a hypergeometric 
construction was proposed for studing arithmetic properties of the values of Dirichlet's 
beta function /3(s) at even positive integers. The construction gives some bonuses [RZ], 
Section 9, for Catalan's constant G = /3(2), such as a second-order Apery-like recursion 
and a permutation group in the sense of G. Rhin and C. Viola [RV]. Here we prove 
expected integrality properties of solutions to the above recursion as well as suggest 
a simpler (also second-order and Apery-like) one for G. We 'enlarge' the permutation 
group of [RZ], Section 9, by showing that the total 120-permutation group of [RV] 
for C(2) can be applied in arithmetic study of Catalan's constant. These considerations 
have computational meanings and do not allow us to prove the (presumed) irrationality 
of G. Finally, we suggest a conjecture yielding the irrationality property of numbers 
(e.g., of Catalan's constant) from existence of suitable second-order difference equations 
(recursions). 

Recently, T. Rivoal and the author [RZ] proved several partial results on the irra- 
tionality of the numbers 



We did not succeed in proving the (expected) irrationality of Catalan's constant 
G = /9(2). However, the general analytic construction in [RZ] allows one to derive 
a certain Apery-like second-order recursion for Catalan's constant; this was done by 
semi-human application of Zeilberger's creative telescoping in [Zu2] and completely 
automatically, thanks to Apery's 'acceleration de la convergence' approach, in [Ze]. 

I would like to thank S. Fischler, T. Rivoal, and J. Sondow for suggestions that 
allowed me to improve the text of the article. 
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1. Hypergeometric series. Recall the Q-linear forms in 1 and G constructed in 
[RZ] and [Zu2]: 



t=0 



(n;=o(i+^'+i))' 



r(3n + 2) r(n + i)^r(?i + 1) 



lim |r„p/" = 



8 r(2?i+ f)3r(2?i + 1) 

3n+l,|n+|, n+^, ^+^5 ?i + 1 

|n + |, 2n + |, 2n + |, 2n + |, 2n + 1 
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lim wl/^ = lim vi/' 



(1) 



Note that Theorem 1 in [Zu2] states the following: 



(2) 



where D^v denotes the least common multiple of the numbers 1, 2, . . . , A'', although 
the better inclusions 



hold for n — 1,2, ... , 1000 by numerical verification (see [Zu2], Section 4). The aim 
of this section is to prove (at least asymptotically, as n — oo, i.e., sufficient for all 
practical purposes) this experimental observation. 

Theorem 1. For n = 0, 1, 2, . . . , we have 



24n+o(n)^^ e Z, 2^^+°(^)L»^^_iV„ G Z. 



(3) 



Rem,ark. As follows from the proof below, the o(n)-term in the inclusions (3) is of 
order log2(2n). 

Proof. We will require Whipple's transform [Ba], Section 4.4, formula (2), 



c, 



^ |a, 1 + a — 6, 1 + a — c, 1 + a — d, 1 + a 



a, 1 + |a, 6, 

r(i + « -rf)r(i + « -e) 

r(l + a)r(l + a-d-e) 
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1 + a — 6 — c, cZ, e 
1 + a — 6, 1 + a — c 



(4) 
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provided that Re(l + a — d — e) > 0, and Bailey's transform [Ba], Section 6.4, for- 
mula (1), 

a, 6, c, d 



4F3{ 1 

k — b, k — c, k — d 

r{k-b)r{k-c)r{k-d) 
r{b) r(c) r{d) r{k-b- c) r{k-b- d) r{k-c- d) 

T{b + t) r(c + t) V{d + t)V{k-a + 2t) 
1 xT{k-b-c-d-t)T(-t) 

J —too 



X 



dt, 



(5) 



27rz J_ioo T(k-a + t) V{k + 2t) 

where the path of integration is parallel to the imaginary axis, except that it is curved, 
if necessary, so that the decreasing sequences of poles of the functions T{k—b—c—d—t) 
and T{—t) lie to the left of the contour, while the increasing sequences of poles of the 
functions V{b + t),V{c + t),V{d + t), QxidV{k - a + 2t) lie to the right. 

Applying transform (4) with a = 3n-|-l, b = c = d = n + \, e = n + 1 and then 
transform (5) with a = 2n + 2, b = n + ^, c = d = n + 1, k = 3n + ^ we obtain 



rr,. = 



(-1)" (2n + 1)! 1 



8 nP 2TTi 

(-1)" (2n + 1)! 1 



r(n + 1 + tfV{n +\ + 2t) V{-tf 
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dt 



T{n+l + t)^T{n+^ + 2t) / tt \ 
r(l + t)2r(3n + f +2t) VsinTTty ^' 



(6) 



Shifting n + 1 +t ^ t and considering the residues at the increasing sequence of poles 
of the integrand in (6) (cf. [Ne], Lemma 2) we arrive at the formula 

oo 



(-1)^+1 (2n + 1)! ^ ^ / T{t)^T{-n - f + 2t) 
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n 



12 



dt \T{-n + t)^T{n + | + 2t) 



-E 



dRnit) 



i/=l 



dt 



t=v 

(7) 



where 



Rn{t) = 



{-ir (2n+i)! (n;=i(^-i))' 



E 

/=0 



A, 



n 



+ 



t-\-l — - — - f-\-l — - — - 



(8) 



The coefficients Ai and ylj, Z = 0, 1, . . . ,n, in the partial- fraction decomposition of 
the function Rnif) can be easily determined by the standard procedure: 



(2n + 1)! 



16 (20!(2n-2i + l)! 

(2n+l)! 



{t-l){t-2)---{t-n) 



{t-l){t-2)---{t-n) 



16 (2Z + 1)! (2n-20! 

Z = 0,l,...,n, 
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hence 

yli-2^"+^eZ and yi;-2^"+^eZ, l = 0,l,...,n, (9) 

by well-known properties of the integer- valued polynomials {t—l){t — 2)---{t — n)/n\. 
Using this decomposition we can continue formula (7) as follows: 

n oo ^ n oo ^ 

rn = i6yAi-y—^-^ + i6yA[-y—^—— 

1=0 IJ.=0 V ' / /^o /Li=0 ^ ^ ' 

m— 1 m—l ^ m' — 1 m' —I ^ 

n I— m — l ^ n ^ — m —1 ^ 

E ^' E E -4! E 

where m = [(n -|- 1)/2J, m' — [n/2\; e = 1 for n even and e = 3 for n odd; e' = 4 — e; 
[ • J denotes the integer part of a number. (For instance, if n is even, we have 



oo n . 

- 

v=l 1=0 ^^^^ 2 4) 

2m 00 ^ 2m oo ^ 

m — l ^ —1 \ 1 00 ^ 



2 



2m / 00 l—m — lx -. 

+ « E -^'(E- E )tz^^ 

l=m+l V=0 /^=0 
2m 00 ^ 

m—l m—l ^ 2m Z— m— 1 ^ 

+ i«E^^E(i^-i« E E 



and we proceed analogously in the three remaining cases.) 

By (8), Rn{t) = 0{t-^) as t ^ oo, hence Er=o^« + Td=o^l = 0- Therefore, 
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formula (10) can be written in the desired form — u^G — f„, where 



Ur 



16(-1)-^A^ = 16(-1)-+i5^A;, 

/=0 1=0 

m—l m—l ^ m' — l m' — l ^ 

/=0 M=l ^ ^ ^ 1=0 n=l ^ 

n I— m — l ^ n l — m' — l 

+ 16 E E 7Z;7T7\2+16 E E 



(11) 



l=m+l n=0 ^ ' > l=m'+l fj.=0 ^ ' > 

(Again, Zeilberger's creative telescoping applied to the sequences u^.v^ in (11) yields 
the Apery- like recursion from [Zu2] for the old sequences w^, in (1); this fact implies 
the coincidence of the two representations (1) and (11) for the numbers and 
in the sequence = u^G — f„.) Formulae (11) for u^^v^ and relations (9) imply 
2^^Uji,2^^Dl^_^Vn G Z. Finally, using (2) and the fact that the order of 2 in Dn is 
Llog2 N\ we arrive at the desired inclusions (3), and the theorem is proved. 

2. A new Apery-like recursion for Catalan's constant. The recursion in [RZ], 
[Zu2] allows one to do fast computation of G with high accuracy. Interpreting the 
solution to the recursion in [RZ], [Zu2] as in (7) prompted us to modify slightly the 
parameters of the above construction. Thus we take the sequence 



= UnG -Vn = - ^ 



dRn{t) 



i/=i 

rn-l 



(12) 



Rn{t) 



2 (n-1)!^ n?:o(2t+i-n-i) 



and apply Zeilberger's algorithm of creative telescoping in order to prove the following 
result. 

Theorem 2. The numbers Un and Vn satisfy the second-order recursion 

{2nf{2n + l)^(20n^ - 20n + 3)wn+i - (3520n^ - 2672n* + 196n^ - 9)w„ 
- {2nf{2n + l)(2n - 3)(20n^ + 20n + 'i)un-i = 0, n = 1, 2, 3, . . . , 

(13) 

with the initial data uq = ^, u\ = 6, and vq = —1, vi = b. In addition, the limit 
relations 



lim \ur,,G - v„\^l'^ 



n— »oo 
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hold and 

24n+o(n) e Z, 2^"+°(")Di^_ii;„ e Z for n= 0,1,2,... . (14) 



(The proof of the inclusions (14) is a word-by- word repetition of what was done in 
Section 1.) 

The polynomials in (13) are polynomials in 2n with integer coefficients; the recur- 
sion (13) looks a little simpler than in [RZ], [Zu2]. As in [Zu2], Theorems 2 and 3, 
the constraint (12) also leads to the continued-fraction expansion 

6G = 5 -I- |-— -J -I- r^-r' -I- r^--' -\ h I TTT^— ' H , 

1^(2) |g(4) |g(6) | g(2n) 

p{n) = (5n2 - 20n + 18) (n - 2)(n - lfn'^{n + lf{n + 2){5n'^ + 20n + 18), 

q{n) = 55n^ - 167n^ + 49n^ - 9, 
and to the multiple Euler-type integral 

«„G-«„ = ^^yj^ ^^^^^ d.d,, 

n= 1,2,3, ... . 

3. A permutation group related to Catalan's constant. Take the parameters 
ho, hi, h2, hs, h^ satisfying the conditions 

hojh^eZ, hi,h2,hs E Z+ ^, (15) 
hj>0 and 1 + ho - hj - hi > for j, I = 1,2,3,4. (16) 

As shown in [RZ], Lemma 2, the quantity 

r(i + ho) r(/i3) r(/i4) r(i + ho-hi- hs) r(i + ho-h2- h^) r(i + ho-h^- h^) 



r(i + ho- hi) r(i + ho- h2) r(i + ho- hs) r(i + ho- h^) 

/lo, 1 + |/io, hi, /i2, hs, h^ 

\ho, 1 + ho - hi,l + ho - h2,l + ho - h^, 1 -|- /iq - /i4 



-1 



(17) 



belongs to the space QG + Q. (In [RZ] , a different ratio of gamma factors multiplies 
the 6^5-series, but one ratio is a rational multiple of the other.) 
By means of the new parameters 



ai = l-\-ho — hi — h2^ a2 = hs, as = h4, 
b2 = 1 + ho - hi, bs = l + ho-h2 
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and thanks to Whipple's transform (4) we can represent the quantity (17) as follows: 

^) 

•1 

dx dy. (18) 



r(a2) r(a3) T{b2 - 02) r(63 - 03) p / oi, 02, as 



r(62)r(63) V ^^2,^3 

1 rl a;a2-l(l _ 2.)62-a2-l^a3-l(i _ ^^bs-as-l 




/o Jo (1 - xyY° 

Finally, take the third 10-element set c: 

Coo = (&2 + ^3) - (oi + 02 + 03) - 1, 

6; — ttj — 1 11 t = 2, 3, 
(hence all c-,7 > —1 by (16)), in order to get that the double integral 

lies in QG + Q. It will be useful to split the set (19) as c = (c', c"), where 
c' = (coo, C21, C22, C33, C31) and c" = (cn, C23, C13, C12, C32) 

will be interpreted as cyclically ordered sets (i.e., cqo follows C31 in c' and cn follows 
C32 in c"). Obviously, each element in c" can be expressed in terms of elements in c', 
and vice versa. Using relations (15) and summarizing what we said above we obtain 
the following result. 

Suppose that 

Coo, C21, C33 e Z + i and C22, C31 G Z (21) 

/or t/ie elements in c' (or, equivalently, C13, C12, 0326^ + ^ arid cn, C23 G Z /or ^/le 
elements in c") and that all elements in c are > —1. Then H{c) e QG + Q. 

Digressing from the demi-integrality of the parameters c, let us note that the hy- 
pergeometric 3F2-representation (18) and the equivalent e^^s-representation (17) lead 
to the following group structure (cf. [Wh] or [Ba], Sections 3.5-3.6). Each permuta- 
tion of the parameters ai, 02, 03 in (18) or of the parameters hi,h2, hs, h4 in (17) gives 
a hypergeometric series of the same kind (but with a different ratio of gamma factors 
before it). For instance, the transposition f) = [hi h^) rearranges the parameters a 
and b as follows: 



ai,a2,a3\ ^ / ^3 - 03, ^2, ^3-^1 

^2, H J V 62 + 63 - ai - 03, 63 



8 



W. ZUDILIN 



and corresponds to Thomae's transformation [Ba], Section 3.2. Hence the group & 
generated by all such permutations appears naturally. An advantage of the superflu- 
ous 10-element set c is the fact that acts on the parameters c quite simply — by 
permutations. As F.J.W. Whipple has shown [Wh], the group is of order 120. 
A possible choice of generators of 6 consists of the transpositions ai = (ai 03), 
CI2 = {0-2 03)5 b = (^2 ^3)5 and the above-cited f) = {hi h^) (see [Zul], Section 6); the 
action of these permutations on the set c reads as follows: 

Cll = (Cii C3i)(ci2 C32)(C13 C33), 02 = (C21 C3i)(c22 C32)(c23 C33), 

(22) 

b = (Ci2 C13)(C22 C23)(C32 C33), i) = (cqo C22)(cil C33)(ci3 C31). 

Theorem 3. Let the quantity H{c) be defined as the double integral in (20), or as 
the ^2-series in (18), or as the oF^-series in (17). Let (25 C ©lo be the c-permutation 
group generated by (22). Suppose that all elements in the set c are > —1. 
Then 

(i) the quantity 

where n(c) = r(coo) r(c2i) r(c22) r(c33) r(c3i), (23) 

is (3 -stable; 

(ii) if the set c is <S>- equivalent to a set satisfying condition (21), we have H{c) e 
QG + Q. 

Proof, (i) The 0-stability of the quantity (23) has to be verified for the permutations 
in the list (22); this is routine using Whipple's transform for verification of the f)- 
stability. 

(ii) In order to deduce the inclusion H{c) e QG + Q from the above claim (i), it 
remains to show that n(crc)/n(c) G Q for a set c satisfying (21) and for all u e 6 or, 
equivalently, for a G {ai, a2, b, [)} (by crc we mean the action of a permutation cr e 
on the set c). This follows easily from the fact that the gamma factors in 

n(c), n(aic), n(a2c), n(bc), n(i)c) 

have exactly three arguments belonging to Z -|- | and two arguments belonging to Z. 

Another (very remarkable) description of the group by means of the double 
integrals (20) and their birational transformations can be found in the work [RV] . 

By [RV], when all elements in c are non-negative integers, one has H{c) G QC(2)-|-Q, 
where C(2) = 7r^/6. Moreover, in this case, D^^Drn^Hic) G ZC(2) -|- Z, where 
"7.1 ^ m2 are the two successive maxima of the set c. This inclusion and the 0- 
stability of the quantity H{c)/Il{c) make it possible to deduce a nice irrationality 
measure for ^(2) (for details, see [RV]). 
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Theorems 1-3 allow us to expect a similar inclusion 

2-'^+o(^)Dm,Dm,Hic) eZG + Z (24) 

if the set c is 0-equivalent to a set satisfying (21); here M is the sum of two integers 
in c' = (coo, C21, C22, C33, C31) and nii ^ 1712 are the two successive maxima of the 
set 2c. Unfortunately, the inclusion (24) is beyond the reach of even the powerful 
group-structure approach to proving irrationality results developed in [RV] (see also 
[Zul]). 

4. Difference equations and irrationality. Since 

lim Dli\ = 

n— >oo 

by the prime number theorem, Theorem 1 (supplemented with equation (1)) or The- 
orem 2 do not yield the irrationality of Catalan's constant. What is the connection 
between irrationality and Apery-like difference equations? We would like to conclude 
this note by pointing out the following expectation. 

A sequence {xn} = {a^nl^o C Q is said to satisfy the geometric condition-^ if the 
least common denominator of the numbers xo,xi, . . . , Xn grows at most geometrically 
as n ^ 00. 

Given a second-order recursion 

Xn+i + a{n)xn + h{n)xn-i = 0, lim a{n) = ao G Q, lim h{n) = bo E Q, (25) 

n— »oo n— »oo 

suppose that the characteristic polynomial + oqA + bo has roots Ai and A2 satis- 
fying < I All < IA2I. Perron's theorem (see, e.g., [Ge], Chapter V, Section 5) then 
guarantees the existence of two linearly independent solutions {x^} and {i/n} such 
that 

lim ^ = Ai, lim ^ = A2. (26) 

Conjecture. In the above notation, suppose that both solutions {xn} and {yn} of 
the recursion (25) are rational and satisfy the geometric condition. Then Ai and A2 
are rational numbers. 

This conjecture is trivially true in the case of constant coeflBcients a{n) = oq and 
b{n) = bo of the recursion (25); we leave this observation as an exercise to the reader. 

In order to show how the irrationality of G follows from the above conjecture, we 
have only to mention that, if G is rational, the solutions {un} ^ind {f„} = {unG — Vn} 
to the recursion (13) are also rational numbers satisfying the geometric condition and 



■"^We should replace the standard term 'G-condition' by the phrase 'geometric condition' since the 
capital letter G is reserved for Catalan's constant here. 
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form Perron's basis, while the roots (ll±5-\/5)/2 = ((l±v^)/2)^ of the characteristic 
polynomial are clearly irrational. 

The geometric condition cannot be removed from hypothesis of the conjecture^. 
Indeed, taking A — (11 + 5v^)/2 and Ai = — 1/A, A2 = A, set 

^n = ^^GQ, yn=LA"JeZ, n = 0,l,2,.... (27) 

Then ~ A^^ and ~ A2 as n ^ 00, hence relations (26) hold. In addition, the 
sequences (27) satisfy the recursion (25) with 



[A-+1J [X^-^P +[X^P'- 
LA"J ^, , [A"J 



n = 0, 1,2, 



[A--1J ^ ^ LA-+1J ' 
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■^It is possible that the conjecture is true if we replace the geometric condition hypothesis by the 
assumption a(n), b{n) 6 Q{n); however this new conjecture would not cover several known cases (for 
instance, the recursion corresponding to Nesterenko's continued fraction for (^(3) in [Ne], Theorem 2). 



